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ABSTRACT 


We consider the acceleration of solar flare protons by cyclotron 
damping of intense Alfven wave turbulence in a magnetic trap. The 
energy diffusion coefficient is computed for a near-isotropic dis- 
tribution of super-Alfvenic protons and a steady-state solution for 
the particle spectrum is found for both transit-time and diffusive 
losses out of the ends of the trap. The acceleration time to a char- 
acteristic energy 20 Mev/nucl can be as short as 10 sec. On the 

basis of phenomenological arguments we infer that the Alfven wave spec- 

_2 

trum has a w frequency dependence and that the correlation time of 

-4 -2 

the turbulence lies in the range 5xl0s<r .^<5x10 s. 

^ corn 


I. INTRODUCTION 


Solar flare protons are believed to be accelerated in complex 
magnetic bipolar configurations in active regions. Evidence points to 
the idea that acceleration occurs in a second stage process [De Jager, 
1969] in contrast to primary energy release of the flare in a first 
stage process which is manifest directly in quasi-thermal heating and 
non-thern’’l 10-500 Kev electrons [Svestka, 1976]. The association of 
Type II radio bursts with proton flares is compatible with the scenario 
that a noisy MHD shock is generated soon {^3 mins) after the initiation 
of the flash phase and accelerates protons by stochastic processes as it 
propagates upwards and away from the flare site [McLean et al., 1971; 

Lin and Hudson, 1976]. The shock is certainly energetic enough and 
current belief is that Fermi processes are at work in conjuction with 
wave-particle scattering to energize protons up to energies > 10 Gev 
[Melrose, 1974; Svestka, 1976]. A large part of our theoretical under- 
standing of acceleration mechanisms has derived from efforts directed 
towards understanding the energization of galactic cosmic rays [Rosen, 
1969]. 

In this paper we will attempt to model the solar proton acceleration 
process with a very simple field-plasma configuration which might be 
described as a "leaky trap". Energization takes place in the corona 
in closed magnetic field structures but allowance is made for particle 
losses. As a first step in this direction we have considered only 
particle escape along the flux tube to the low solar atmosphere. The 
acceleration mechanism is cyclotron damping by super-Alfvenic protons 
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of intense Alfven waves propagating in the trap. This process is 
intrinsically stochastic and no adiabatic, Fermi-type mechanism is 
invoked for the acceleration. 

Davis [1956] gave the earliest description of the effect that a 
turbulent state of MHD waves produces a diffusion in pitch angles 
towards isotropy. The development of quasi-linear theory [see Kennel 
and Engleniann, 1966, and their refs.] provided a useful analytic tool 
with which one could address the problem of diffusion in phase-space 
resulting from turbulent wave-particle interactions. Sturrock [1966], 
using a Fokker-Planck formalism, achieved similar results and thus 
fortified our confidence in these methods. 

Melrose [1974] demonstrated that the combined effects of stochasti 
pitch angle scattering by Alfven waves and Fermi acceleration could 
account for the energization of solar flare protons. The source of the 
Alfven turbulence was considered to be the non-thermal protons them- 
selves which, if sufficiently anisotropic in pitch angle distribution, 
would drive an electromagnetic Alfven wave instability. Since the free 
energy for the waves is derived from the kinetic energy of the non- 
thermal protons in this picture, another, independent, free energy 
source must energize the protons (e.g., a Fermi process). In our model 
the energy flow is opposite to that of Melrose. We assume that Alfven 
turbulence is generated ad hoc by an exciting agency that is flare 
related. Using this approach we may eliminate the Fermi mechanism and 
consider the acceleration possible from just short wavelength, kv = 


Alfven waves which not only scatter in pitch angle but also energize 
through rollisionless damping [Schatzman, 1967’, Sturrock, 1974]. 

The model can be applied to the system described by, for instance, 
McLean et al. [1971] in which a localized shock generates Alfven waves 
propagating throughout a flux tube. With some modifications, the model 
may also be applied to a system like that of Sturrock [1974] in which 
protons are accelerated on open field lines by a noisy MHD shock. 
However, recent observations have suggested that some protons may 

be accelerated as early as the impulsive phase. We have, therefore, 
concentrated on the consequences of a leaky trap without pinning down 
the source of the turbulence or when and how it is generated. If 
impulsive phase proton acceleration is confirmed, then we hope our 
model may provide insight into the primary energy release mechanism 
itself. 

In Section II we outline the ingredients of the model and inquire 
into the possibility of proton acceleration during primary energy 
release. In Section III we derive the energy diffusion coefficient, 

for the cyclotron damping mechanism for a near-isotropic distribu- 
tion of super-Alfvenic ions. In Section IV we apply this result to 
the leaky trap with either ballistic transit-time losses or diffusive 
losses out the ends of the trap. In the final section we summarize 
our results and comment on the possibility of flare-generated MHD 
standing waves in the trap. 
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II. MODEL 


A. Basic Assumptions 

The basic magnetic field configuration we shall consider is that of a 
closed loop, bipolar arch structure. We presume that some exciting 
agency, that is flare related, generates intense MHD turbulence which 
interacts with particles at coronal heights. The closed field lines 
serve as a magnetic trap for the non-thermal particles which are stochas- 
tically accelerated by the turbulence but are prevented from escaping by 
mirror forces arising from the convergence of the mean field towards 
the photosphere. The turbulence is assumed to consist of short wavelength 
perturbations, X << L, much smaller than any macroscopic scale of the 
system; this enables us to use homogeneous quasi-1 inear theory to gauge 
the acceleration efficiency. The dominant interaction with the waves is 
considered to be with the wavevector component parallel to [i.e., 
ki= o] and we, therefore, model the turbulence as an ensemble of parallel- 
propagating Alfven waves in the frequency range < w << where 
is the proton gyrofrequency and is a low frequency cutoff. 

The question of particle phase-space losses is an important one and 
different assumptions concerning what loss mechanism prevails will alter 
the picture significantly. We shall consider proton energies sufficiently 
large that they may be considered collisionless with respect to binary 
coulomb interactions [Braginskii, 1965] during the acceleration period. 
However, if the protons remain in dense regions sufficievitly long, this 
deenergization process must be accounted for. We may expect that some 
cross-field diffusion of particles out of the trap is likely to occur 
during the acceleration period. However, we have chosen to consider what 
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may be the fastest loss process; particle escape along the field to the 
chromosphere and photosphere, A consequence of our assumption of intense 
Alfven turbulence (with a predominant magnetic component) is that particles 
will be scattered very efficiently in pitch angle and can either stream 
out or diffuse out of the acceleration region along the flux tube to the 
lower atmosphere. For either case we have scaled the losstime to the 
particle ballistic travel time and the diffusion time along the length of 
the flux tube. In the first case (transit-time losses) our model is anal- 
ogous to the process described by Kennel [1969] who considered stochastic 
acceleration of electrons in the earth's magnetosphere. In the second 
case (diffusive losses) the model incorporates diffusive spatial transport 
in the presence of elastic wave collisions [Jokipii, 1971] analogous to 
the propagation effect observed in interplanetary space. 

B. Applications 

The association of Type II radio bursts with proton events observed 
at earth suggests that a flare-initiated shock propagates upwards and 
accelerates protons. If the shock is the exciting agent of MHD waves, 
then our model may be used to compute the acceleration of protons in a 
configuration like that proposed by McLean et al . [1971]. It must be 
determined, for instance, if the acceleration takes place ahead of or 
behind the shock, and the question of how the protons obtain access to 
interplanetary field lines must be addressed. If the field lines of the 
trap "open up" [Kopp and Pneuman, 1976], then particles will be injected 
ball istically into the solar wind; but if the field lines remain closed, 
but distended, cross-field diffusion effects must occur to exhaust the 
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trap. Chronologically, this process is an “after effect" of the flare; 
that is, subsequent to or during the primary energy release the shock 
is formed and proton acceleration begins at times ^ 3 mins after the 
initiation of the flash phase. 

However, y-^’ay observations [Chupp et al . , 1973] made by the OSO-7 
satellite have indicated that proton acceleration may begin at times coin- 
cident with the flash phase onset. Figure 1 is a plot of the Aug. 4, 1972 
2.2 Mev Y-i^^y line together with hard x-ray profiles from Van Beck et al . 
[1973]. We note that while the maximum of the y-^ay i^lux occurs after 
the hard x-ray maximum a not insignificant flux is measured very soon after 
onset [Svestka, 1976]. In fact, the flux at 0625 UT is 1/3 of the maximum. 
Since there must be a finite acceleration time for the protons to reach 
these energies in sufficient number. Figure 1 indicates that, within the 
temporal resolution of the instrument and ignoring any nuclear reaction 
time delay, a generous upper limit for the acceleration time a- 3 mins 
(assuming that energization begins at 0622). The acceleration time could 
be much less (lOs < < 200s) and still be consistent with the non- 

thermal 2.2 Mev emission measure reaching a maximum later as more flux 
tubes become involved and the acceleration region increases. The measure- 
ment at 0622, if real, gives support for a short ^-lOs acceleration time. 
Ramaty et al . [1*^75] interpret a break in the hard x-ray spectrum of 
Aug. 4, 1972 at about 500 Kev towards a harder profile at higher energies 
as evidence for electrons undergoing second stage acceleration. 
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Since we might expect intense Alfven wave turbulence generated in 
primary energy release mechanisms like, for instance, magnetic field 
tearing instabilities [Spicer, 1976], we have, therefore, considered 
a parameter regime for the leaky trap that might be applied to plasma- 
field conditions in the vicinity of x-ray kernals and their coronal 
loops at the time of the impulsive phase. The model, then, can allow 
us tG infer something about the primary energy release process itself 
rather than post-explosion conditions in MHD shocks. Future observations 
should define clearly the role of magnetic turbulence, especially with 
regard to primary energy release, and can test our supposition that 
proton acceleration can occur as early as the impulsive phase. 
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III. DIFFUSION COEFFICIENT 


A. Quasi -I inear Theory 

We consider a homogeneous plasma in a magnetic field directed 
along the z axis and assume that the turbulence, composed of parallel-* 
propagating Alfven waves {k^ = o) with frequencies in the range 

“o ' “ = "cp- 

can be described by a spectral density 


< ^1/Sir > = 2 / du/^(b)). 

O 


( 2 ) 


= B^/(4TrMnQ)'^ is the Alfven velocity, = eB^/Mc is the proton 
gyrofrequency, and is the perturbation magnetic field related to the 
induced electric field by 

|6, 1 = ''i l«, I • (3) 

The factor of two in equation (2) is to discriminate between waves travel- 
ling parallel and antiparallel to 

It is well known that the temporal evolution of a distribution of 
particles, say protons, in response to resonant scattering by Alfven waves, 
can be described by a quasi -linear relation [Kennel and Engleman, 1966; 

Hall and Sturrock, 1967; Davidson, 1972] which, for ki. = o, reduces to 
in the non-relativistic regime 






Vx 



(4) 


where 


J = 


4ir^e^V^ 

M^c^ 



3f 

3Vx 


+ 



(5) 
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In the frequency' range considered, the group velocity is the Alfven speed, 
[Vgj = V^; and v^ are unit vectors in velocity space (cylindrical 
coordinates}, azimuthal symmetry about being assumed. The effect of 
the finite perturbation electric field (3) is represented in the terms 
consisting of the integer 1. The unique relationship between the proton 
and the Alfven wave in this interaction is expressed in the cyclotron 
resonance condition which is given as 

t. - kv^ = ±fi^p . (6) 

the upper (lower) sign taken for the L (R) mode of the Alfven wave. For 
k > 0 the terms L and R also give the polarization of the wave. Equation 
(6) indicates that an effective "collision" occurs when the doppler- 
shifted wave frequency is the proton's gyrofrequency and the sense of 
rotation is the same. 

For the range of parameters we shall consider, « car'd the effect 
of the predominantly magnetic wave is to scatter the particles primarily 
in pitch angle. It is convenient, therefore, to express (5) in spherical 
coordinates. If cose = v^/v we have, using (6), 
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( 7 ) 


Inspection of (6) indicates that if oj << Si , then |v | >> V» if 

Cp Z A 


resonant scattering is to occur. Protons with |v^| < follow the wave 


adiabatical ly conserving their first invariant \i = during the 
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oscniation, whereas super-Alfvenic protons can satisfy (6) and are 
non-adiabatical ly scattered, the interaction being sufficiently rapid 
(t -V so as to violate their first invariant. Equation (6) then 

cp 

permits a useful perturbation expansion [Melrose, 1974] 

_ +1/VZ _ 7 


k/si 




“z ""z 

m, — 


( 8 ) 


for particles with velocities |v^| » V^. Retention of the zero order 
term in (8) provides strictly magnetic effects, while the first order 
term incorporates effects of a finite electric field due to oi o. The 
expansion of (7) to zero order yields 

M^c2 I cose I ® 

Any reference to frequency here is incidental since if (2) is given 
alternately as 




( 10 ) 


we have the relation 


=A(k) (11) 

and thus (9) can also represent non-adiabatic pitch angle scattering from 
stationary magnetic inhomogeneities descrioed byfi(k) [Jokipii, 1971]. 

On the shortest time scale, the turbulence acts to isotropize the 
distribution and on a longer time scale O(vVV^) the particles are ener- 
gized through stochastic accelerations from . This results in a diffu- 
( 2 ) 

Sion current ' which we may conveniently write using the approximation 
that the distribution is nearly isotropic from pitch angle scattering and 

ri f 

^ = 0 . Retaining the first order term in (8) provides the second order 
term for = 0(V2/v2) 
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j(2) . i’!!!!)! 


V.sin^e , . 


M^c^- V I cose 

and the distribution evolves as 
4Tr2e^V| 




v|i] 

Icosel v^ 9v L J 
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B. Wave Spectrum 

To model the turbulence we assume the spectral density can be 
represented by a power law in frequency 

. (14) 

o(x) is the Heaviside step function which introduces a low frequency 
cutoff at ojp. It is also implicit in (14) that frequencies do not extend 
as high as the proton gyrofrequency. The resonance condition (6) gives 
the inverse wave-particle relation 




(15) 


V I cose I 

The subsequent steps are straightforward: the spectral density is 

normalized to the perturbation field energy density in (2) to determine 
equations (14) and (15) are inserted into (13) and the equation is 
averaged over pitch angles, and the change of variable to energy per nucl 
is made. 

Generalizing this procedure for an arbitrary ion of ionic charge Ze 
and mass AM, the temporal evolution of the differential number density for 
arty ionic species can be written as 


at 


N(E,t) = D ^ [ E 


n+3/^ a 




aE 


/E 


(16) 
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where the total number density of the ionic species is given as 


and 


N = TdE N(E) 

0 

D = (A/Z)^*'’^ D 


P ’ 


P B2/8ir (2n+l)(2n+3) e" 


(17) 

(18) 

(19) 


mp 


B, 


2 2 “0 
The variable E = ^jMv = energy/nucl of the species. E^ = %MV^ = g— — 

is the Alfven energy, M = proton mass, and the index n is related 

to the wave spectral index by 

n = > 5 ( 6 - 1 ). (20) 

The normalization of A(w) introduced the low frequency cutoff 

into the scaling of the diffusion coefficient and we have written it in 

energy units as 


£ = UMV? f 
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( 21 ) 


Ej^p is related to the maximum energy a proton can achieve in the following 
way: if is non-relativistic in value, then E^g^.p ^ ^mp’ W 

is relativistic, one must use the relativistic form of the resonance 
condition instead of (6) to determine E . Thus, the maximum energy 

Inc] A 5 p 

a proton can be accelerated to is 

E, 

I 

E 


E„_ « IGev 
mp 


■max,p 


mp 

V. 

Mc 2 ^_^ : E„_ » IGev 


( 22 ) 


c 


mp 


Again, so long as we are considering the non-relativistic regime of 
energies in (16), (21) is used by definition for all values of 

It is of interest to ask what minimum energy (injection energy) of 
electrons is necessary for them to be accelerated by Alfven waves. The 



resonance condition for electrons with gyrofrequency = eB^/mc is 
given as 


- kv^ = + s^ce/^ 


( 23 ) 


which indicates that if t*' < n , we have [Sturrock, 1974] 


^inj,p '~l\ 


= E„ - 


m “"A 


non-rel 


(24) 


^inj.e 


C Mc^ ~ ; ultra-rel 

These restrictions apply for MHO Alfven waves where < n . Extension 

cp 

of the wave spectrum to higher frequencies to include whistler waves 
(R mode) or Langmuir waves lowers the electron energy threshold [Melrose, 
1974], but is better addressed as a separate problem since the wave 
properties are quite different and, phenomenologically, acceleration by 
higher frequency waves may be related to first stage acceleration 
processes [Hoyng, 1977]. 

If, instead of energy, one considers a diffusion in the non-relativis- 
tic regime of rigidity, R = AMvc/Ze, we have alternatively 


where 


it N(R.t) = K ^ [ r2-"^2 ^ 3 

. Ai h/8iy ^ 

(2n-')(2n.3) R^J 


(25) 


(26) 


C. Physical Interpretation 

An interpretation of the physical process leading to equations (4), 
(6), and (14) is given in Figure 2. We have decomposed the turbulence 
into an ensemble of elemental scattering centers of length x < 

This ensemble is equivalent to the actual situation in that the idealized 
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scatterers have the same autocorrelation function as the actual system, 

p{t) = ] 

p(x) =./dkfl(k)e’'<V[i"dk4(k) ] . 

From the power spectrum given in (14) we identify (6 > 1) 

^corr ~ ^ ^corr " ^corr ’ 

0 

and d(K) [orfi( (u)] is a measure of the differential number density of 
elemental scattering centers of size k"^ 's x < x . A super-Alfvenic 

V w I r 

proton satisfying the resonance condition (6) can be effectively scattered 
as drawn in Figure 2. Particle A sees a component of parallel to v 
throughout the interaction and is accelerated, while particle B is 180° 
out of phase and is decelerated. A collection of random (cyclotron) 
phased particles will experience, however, a dispersion in energies. 

The acceleration process described here cannot be interpreted as a 
stochastic Fermic type. The acceleration mechanism described by Fermi 
[1949, 1954] involves an adiabatic process such that the particles' first 

mw2 ^ 

and second invariants, y = -^/2B, 0 =©Vj|ds, are conserved. Pitch 

angle scattering is an ancillary process invoked to counterbalance the 
systematic decrease in the particle's pitch angle as it is accelerated 
[Davis, 1956]. We have described a process where oscillations are suffi- 
ciently rapid that they non-adiabatical ly violate the first invariant. 
Energization (or denergization) is concomitant with pitch angle scattering 
when ui f 0 and such a process is best interpreted as cyclotron damping 
of Alfven waves [Stix, 1962]. 


(27) 

(28) 

(29) 
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IV. CHARACTERISTIC SCALINGS 

In order to apply the result of equation (16) to the leaky trap 
model, we have considered three problems: a transient solution and a 

steady-state solution with transit-time and diffusive losses. 

A. Initial Value Problem 

Two facts which immediately emerge from injection of (16) is that 

1 „ 

the energization time scales as E and that the solution asymptotes 
to a "steady state" N(E) E"*^. The energization time indicates that 
n = 1 distinguishes two, quite different, diffusion regimes which favor 
acceleration of either high or low energy particles. We shall argue in 
this paper that the regime which favors low energy particle acceleration 
is more appropriate for solar proton flares and we shall consider the 
parameter regime 

1<6<3 0<n<l. (30) 

This statement is equivalent to the expectation that if the power spec- 
trum in (14) is sufficiently flat, protons can be accelerated over a 
large energy range in the shortest time for a given perturbation field 
energy density. Increasing 6 to larger values puts all the energy in 
the low frequencies and biases the acceleration in favor of very high 
energy particles [Equation (15)]. 

The homogeneous solution of (16) with the boundary conditions at an 
injection energy E^^ that 

N(Eo,t) = Nq , N(E>EQ,t=0) = 0, N(E ^ «) -v 0 
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is given as 


N(E 


/ la. 

s 


K.feE” /i") 


ds gS ( E )P Iji 

^ ^ ^ K (,e:: ’ 


( 31 ) 


- 1 ® 


with 


P = k{}-2n) 
w = ^(1-n) > 0 


K = 


ri 


^JVot' 


_ 2n+l 


(32) 


If we consider times long enough such thatsE^ « 1, the small- 
argument form of K^(^) can be used in the denominator and the integral 


is evaluated as 


N(E,t) = Nq 


n r{y, E2-W) 
r(y) 


(33) 


where r(a,x) is the incomplete gamma function [Abromowitz and Stegun, 1970]. 
This function exhibits two different behaviors depending on whether 


X = Js5 2 e2W > 1 , 

We may approximate the solution as 


(34) 


N(E,t) 



N. 


Eo -,n 


0 

E 


n (i« 52 e 2 W)»^’'’ 




oW 


r(y) 


X « 1 


X » 1 


(35) 


The solution changes from a power low at low energies to an exponential 
behavior (w>0) at the diffusion front characterized by a break energy. 
Eg, which depends on time and is defined as x = 1 leading to 


1 


Eg(t) = [(l-n)2 Dt] 


(36) 
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In a subsequent section we shall argue for the feasibility of the 
value n = Jj. Anticipating this discussion and defining eg = we may 

scale the break energy as 


r - n pi f ^0 Y f f 1 6xl0^cm"3 V 

Eg(n-^,t) - 0.61 ^ ^ J ( — J J 


t^ Mev/nucl . 


Time is in units of seconds and our choice of and n^ = background 
density corresponds to an Alfven energy of Ey^ * 41.4 Kev. Note also 
[Equation 18] that for n = 4. reference to Z/A disappears and (37) is valid 


for any ionic species. The value 10 Gev for E^^p corresponds to p = 
4.3 Gev. 


i 


B. Steady-State Solution 

Assuming that proton acceleration can occur for a period longer than 
10 sec, (37) admits the possibility of very efficient proton acceleration. 
Eg being a measure of the "hardness" of the proton spectrum. A realistic 
model should allow for particle losses and phenomenologically we may write 

] - 0 (38) 

to describe the time asymptotic evolution of the transient solution (33). 
The steady-state solution of (38) should give a meaningful result provided 
that the fastest losstime 

T (39) 


3N _ r, r cf'+3/2 j_ f X 

3t ■ ^ 9E L t 3E i ^ 


I 
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can be ascertained. The solution for q = ^sCl-n-e) > 0 is 

,p kJxe'I) 


N(E) = N„ (f 

^ n J 


0 ^ K ,(XE^) 


( 40 ) 


where 


p = Ml-2n) 


X = 


V’ O' 


(l-n-0) 




q = i5(l-n-|3) > 0 V = - 20-n-B) ‘ ^ 

If XEq « 1 we may use small argument form for the denominator and 
(40) exhibits the following behavior 

: XE^ « 1 


N, (E^/E)" 


N(E) 


N -?£- 
° r(v) 


%XE, 


V-% ' 


0 

E 


!.(n-e) ^-xe'i . Jjq , _ 


(42) 


The steady-state break energy, defined as xE*^ =1, is 

1 


Eg = [ ls(l-n-S)2 0 t„ ] . 


(43) 


C. Transit-Time Losses 

In a magnetic field minor geometry, particles with pitch angles 

less than the loss cone value 

sin^G = 

0 0 max 


(44) 


can penetrate to a field strength of B_„ where they may be considered as 
lost [through significant denergization processes in the dense chromosphere 
and photosphere]. This angle corresponds to a loss cone "volume" of 
4-!r[l-cose^] steradians. If pitch angle scattering occurs rapidly enough 
to maintain a significant number of particles in the loss cone, then the 
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losstime is estimated to be the transit time along the flux tube with a 
weighting factor of the fractional volume of the loss cone [the strong 
pitch angle diffusion regime of Kennel, 1969], i.e., 

T = (M/2)^ L(l-cose^j)"^ E“^ . (45) 


Transit-time losses lead to a value of 3 = in (39). We note that 
if n = hi then the solution (40) behaves at large energies as N(n=-s=-'i) 
Observations made at lAU have indicated that over a large 
range of energies, proton spectra can be fit by an exponential in rigidity, 
the fit improving at high energies (^200MV = 21Mev) [Frcier and Webber, 

1963; Van Hollebeke et al . , 1975]. An exponential in rigidity is also 
consistent with models based on y-ray observations [Ramaty and Lingenfelter, 
1972]. On this basis we argue for the feasibility of the value r\ = h 
corresponding to a wave spectral index of 6 = 2. 

If we take a coronal field strength of lOOG and a photospheric field 
of lOOOG, the loss cone angle = 18°. Using (45) in (43) gives a steady- 
state break energy for transit-time losses of 


Eg (n=-0=}5) = 11 


B 


0 


1006 
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^B 1 

1 f L 1 


f lOGev 1 

1% 

10‘2 i 

1 lO^cm ^ 


^ ^mp ^ 



6xl0^cm"^ 1 


f — 1 Mev/nucl . (46) 

^ o' max ^ 

Although (46) gives a reasonable value for the turnover, the scaling 
may not be correct. In the absence of cross-field losses, (45) without 
the loss cone weighting factor must be the shortest conceivable time to 
deplete the trap. But, for the parameters chosen, the mean free path for 
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pitch angle scattering is much less than the length of the tube. Under 
these circumstances, a spatial diffusion along the field to the lower 
atmosphere may be more appropriate than losses scaled to the particle 
ballistic transit time. 


D. Diffusive Losses 

If the phase-space distribution function varies along the z direction 
and we consider only the effects of pitch angle scattering by waves, then 
(4) in zero order can be written as 


If * V, ^ f (Z,7,t) = ^ t 0. 


8f 


] 


(47) 


9n nn 9n 

and n = cose. Jokipii [1971] has shown that in the approximation that f 
is nearly isotropic, (47) reduces to a diffusion equation for the differ- 
ential number density 

(48) 


iN(z.E.t) = K,, 


where 


■'ll ^ y 

and from (9) we have for protons 
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Substituting the spectral density (14) in the above and generalizing for 
an arbitrary ion yields the result 


(50) 
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(51) 
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and the parallel spatial diffusion coefficient is 
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( 52 ) 
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If we consider ions diffusing down both legs of a trap of length L and 

I ^ 

estimate from (48) the losstime as — , we have 
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mp 


If n = then - -h and the steady state solution will behave in a 
manner similar to the results of the previous section for transit time 
losses. Equations (53) and (43) give a turnover for diffusive losses of 


Eg (n=is) = 17.2 
Mev/nuci . 
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V. SUMMARY AND DISCUSSION 

We have considered some of the consequences of a turbulent spectrum 
of intense Alfven waves propagating along the magnetic field of coronal 
flux tubes. The model describes a steady-state situation where a popula- 
tion of ions is maintained at a threshold velocity of the Alfven speed 
and is scattered in pitch angles by the turbulence towards isotropy 
and energized concomitantly through collisionless damping. The accelera- 
tion time to MOMev energies can be as short as -x-lOs and the process itself 
is limited only by how well the trap can contain the accelerated particles. 
The measure of the acceleration efficiency has been gauged by a break 
energy of the particle spectrum which characterizes a turnover from a 
relatively flat partical spectrum to a rapidly decreasing exponential at 
high energies. This break energy is approximately the value at which the 
acceleration time equals the loss time. 

The primary loss mechanisms which we have considered are transit- 
time losses and diffusive losses out of the ends of the flux tube. For 
the range of parameters such that the mean-free-path for pitch angle 
scattering is small relative to the length of the tube, the diffusive 
picture is probably more correct. The question of other losses is impor- 
tant and further work is needed to clarify the role of cross-field diffu- 
sive effects, especially with regard to how particles obtain access to 
interplanetary field lines. However, if the acceleration region consists 
of open field lines [Sturrock, 1974] this problem disappear". 

One important conclusion we have drawn is that for a given r.m.s. 
magnetic perturbation B-| the wave power spectrum in (14) must be relatively 
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flat (30) if efficient acceleration is to occur. For large values of 
the wa'e spectral index 6 the diffusion process favors the acceleration 
of the very high energy particles which are scattered by the most ener- 
getic oscillations near the correlation length of the turbulence; in 
this case the energy diffusion coefficient scales as a strong power of 
(E^/Ej^p)^ [Newman, 1975] indicating that the acceleration of low energy 
particles near the injection energy is inefficient and holding up 

xJ 

the works. On the basis, of this model proton flares are associated with 
both large magnetic field perturbations and flat power spectra (30). 

This result is consistent with the findings of La Combe [1977] who con- 
sidered particle acceleration in radio galaxies by Alfven waves. We may 
also extend the result (19) to include the value & - l(n=0) for which 
the substitution is made in both (19) and (51) 


2n 


1 

(2n+l)(2n+3) 3 AT 

1 n max o 

In particular, we have shown that the value 6=2 leads to an expo- 
nential particle spectrum at high energies that behaves as N E'"^ 

On the basis of observational indications of this behavior we have assumed 
6 = 2 as a likely value. For this case, we may also rewrite the results 
(52) and (54) in a simple form in terms of the low frequency cutoff Uq 
and we offer the following speculations on second stage acceleration in 
conclusion; 


(55) 


The choice of scaling the diffusion coefficients in terms of the quantity 
E^p (?1 ) throughout this paper was motivated by the possibility of deter- 
mining the correlation time, from observations of the highest 
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recorded proton energies. The value taken of E^p = lOGev corresponds 

^max,p “ (21) we infer a correlation time and 

distance of 


corr 


corr 




sec 


( 56 ) 


It is possible, however, that the correlation time is much longer than 
(56) and the much higher proton energies possible are not recorded because 
of the rapid decrease in the spectrum (40) at very high energies. Treating 
uJq then as a free parameter and taking n = Jg, we may rewrite in summary 
the results for spatial diffusion as [e„ = 

b 10 


and 


‘f “o 
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(57) 


(58) 


Defining a mean-free-path such that K, , = i v \ we have for spatial 


diffusion 




o^B 


(59) 


The energy diffusion coefficient (19) becomes 

DOl-y = i Ej E^ 

leading to the temporal break energy (36) of 

EB(n=!s,t) = 


_TT_ 

16 


1 ('/^corr)^ 


(60) 


(61) 
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and a steady-state break energy for diffusive losses of 


■ 2048 “^0 


If we estimate from (16) an energy diffusion rate 
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( 62 ) 

(63) 


we get a result which differs from (61) by a factor of 4 but is equal to 
Melrose's [1974] result [see his equation 24 and the discussion following]. 
Equation (62) may be scaled as 


Eg(n=3i) = 18 


0.05s 12 [ Ib 12 f L 
^corr J t 1 J i lO^cm 


2 

Mev/nucl , 


(64) 


and it is clear that longer correlation times than (56) require a larger 
Eg in proportion to achieve a break energy in the range of 20Mev/nucl . 

In fact if we consider that the value (64) for the break energy is typical 
we may place the correlation time in the range 

5xl0'‘*s < < 5xl0"2s , (65) 


where the smaller limit is determined by the parameters given by (56) such 
that > 4.3Gev and the larger limit is determined by (64) such that 

Eg < 1 . 

The limit Eg = 1 offers an interesting interpretation in terms of low 
frequency MHO standing waves of the trap. If the diffusive scale length 
is L/2 and the wavenumbers of the magnetic perturbations take the harmonic 
values k = Z-nn/!^ where ^ < L, the break energy for diffusive losses (62) 
can be written as 

( X "corr h 

where k^^ = is the dominant wavenumber excited. If we require an 
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energization of about Eg = lO^ E^, and suppose that, say, L = 3j(, then 
the limit = 1 is consistent with a low harmonic, = 6, being 

excited during the acceleration. The fundamental frequency, = 
gives a greatest upper bound to the correlation time of 


T < T 

corr c 


f lOOG 

f "o 1 

!s f i ] 

1 «o J 

' 6xl0^cm"3 ■ 
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( 67 ) 


We have in this paper considered the possibility that second stage 
acceleration may occur during an/or immediately following the impulsive 
phase of proton flares. If primary energy release is associated with 
large amplitude Alfven waves self-consistently generated, then (66) admits 
the possibility of low harmonic standing waves playing a role in the basic 
flare instability. Phenomenologically, we may speculate that the correla- 
tion time (65) is related to the pulse width of hard x-ray bursts. 
Presently, observations show a fine structure to the bursts consisting 
of -V Isec spikes [Svestka, 1976]. The Goddard x-ray spectrometer aboard 
the upcoming Solar Maximum Mission (SMM) satellite will have a temporal 
resolution lO'^ s < x < 10”i s from 20-300 Kev energies. If, indeed, 
there is a relation between the Alfven wave correlation time and the hard 
x-ray pulse widths, we believe that the chances are good that this instru- 
ment will detect fine structure in the range of (65). The University of 
New Hampshire/Max Planck Institute gamma ray experiment also aboard SMM, 
with a temporal resolution down to 64 ms, should settle the question of 
how soon are >10 Mev protons accelerated. 
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FIGURE CAPTIONS 


Figure 1. The temporal evolution of the hard x-ray and 2.2 Mev 7 -ray 
bursts for the August 4, 1972 flare (from Lin and Hudson, 
1976). 

Figure 2. Resonant interaction of super-Alfvenic protons with a linear 
polarized elemental scattering center of length x < 

The solid line portions of the trajectories indicate posi- 
tions on this side of the page. Energization or deenergiza- 
tion depends on the cyclotron phase of the test particle. 

An average over a random- phased, isotropic distribution of 
particles leads to cyclotron damping of the oscillation. 
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